Low-lying excitations of a trapped rotating Bose-Einstein condensate 
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We investigate the low-lying excitations of a weakly-interacting, harmonically-trapped Bose-Einstein 
condensed gas under rotation, in the limit where the angular mometum L of the system is much 
less than the number of the atoms A'^ in the trap. We show that in the asymptotic limit N oo 
the excitation energy, measured from the energy of the lowest state, is given by 21Nz{Nz — l)uo/68, 
where is the number of octupole excitations and va is the unit of the interaction energy. 

PACS numbers: 03.75.Fi, 05.30. Jp, 67.40.Db, 67.40.Vs 



The behaviour of a trapped Bose-Einstein condensate 
under rotation has attracted much attention in the recent 
years. References ^]-|l3| have dealt with this problem 
theoretically both in the Thomas-Fermi limit of strong 
interactions, as well as in the limit of weak interactions. 
Borrowed from the field of nuclear physics, the termi- 
nology "yrast state" refers to the state with the lowest 
energy for a given angular momentum L that the system 
has. Of equal importance to the yrast state are the low- 
lying excitations, which determine the thermodynamic 
behaviour of the system, and they are also crucial for the 
stability of the yrast states under external perturbations. 

In the present study we consider the hamiltonian H = 
Ho + V. Here 
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includes the kinetic energy of the particles and the poten- 
tial energy due to the trapping potential, where M is the 
atom mass, and lo is the frequency of the harmonic con- 
fining potential, which we assume to be isotropic. Also 
the interaction between the particles V is assumed to be 
short-ranged. 
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where C/q = Anxh^ajM is the matrix element for atom- 
atom collisions, with a being the scattering length. We 
consider the limit of weak interactions, where the typical 
interaction energy is much less than the typical oscillator 
quantum of energy. 



nUn <C 'hio. 
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with n being the density of atoms. The above condition 
allows us to work in the subspace of single-particle states 
with no radial excitations. 
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Here p, z, and are cylindrical polar coordinates, and 
o-o = {fi/MujY/'^ is the oscillator length. These states 
are degenerate in the absence of interactions (while states 
with radial excitations lie higher by an energy of order 
TiLj). The whole problem thus reduces to incorporating 
the interactions between the atoms, which lift the degen- 
eracy of the states. 

One of us has determined in Ref. the yrast line of a 
weakly-interacting Bose-Einstein condensate for an effec- 
tive attraction between the atoms. In the same reference 
the yrast line has also been determined for the case of 
an effective repulsion between the atoms and for L ^ N. 
In Ref. we have given a more detailed description of 
the yrast line within the mean-field Gross-Pitaevskii ap- 
proximation for a wide range of values of the ratio L /N. 

If one studies the limit L <^ N, a,s shown in Ref. [||, 
in a state with a 2^-pole excitation. 
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where m is the state with angular momentum mh, the 
excitation (interaction) energy is 



eA - {X\V\X) - {0\V\0) - ~Nvo 1 
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0) ) is the ground state, and 



where |0) — {{m 
vo — Uo/qq. The above equation implies that the high- 
est gain in the interaction energy per unit of angular 
momentum comes from the A = 2 or 3 excitations, and 
therefore the quadrupole and octupole excitations are ex- 
pected to carry the angular momentum for L <^ N. As a 
result, in this limit there is a quasi-degeneracy between 
the low-lying states, as we discuss below in detail. By 
saying quasi-degenerate states, we mean that the energy 
separation between them is of order vo, and not of order 
Nvo. 

Bertsch and Papenbrock Q have examined the ground 
state of a weakly-interacting Bose-Einstein condensate 
under rotation numerically by diagonalizing the hamilto- 
nian H in the subspace of degenerate states (^ for a given 
angular momentum. In Ref. Nakajima and Ueda 
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have performed similar numerical calculations in the limit 
where the angular momentum per particle L/N is much 
less than 1, and have found that the quasi-degenerate 
states which we mentioned in the previous paragraph lie 
above the yrast by an energy, which in the asymptotic 
limit ^ oo is given by 1.59A'3(iV3 — l)i;o/4, where 
is the number of octupole excitations. 

In this study we give an analytical derivation of this 
result with use of a diagrammatic perturbation-theory 
approach. The starting point of our analysis is the fact 
that the quadrupolc and octupole A = 2, 3 excitations are 
dominant for L <^ N |^Jl^ , and we therefore assume that 
in this limit the angular momentum is carried by A = 2 
and A = 3 excitations only. In addition, the condensate 
is dominated by atoms which do not have any angular 
momentum. Our approach therefore consists of consid- 
ering a condensate with A'o atoms in the state m — 0, N2 
atoms in the state with m — 2, and N3 in the state with 
m = 3, and then treating the other states perturbatively 
by keeping the appropriate diagrams and using pertur- 
bation theory to get the correction to the energy. The 
"bare" interaction energy with particles in the m = 0, 2 
and 3 states is given by (see Fig. 1 and Table I) 

f = Ino{No - 1) + ^N2{N2 - 1) + - 1) + 

i lb 32 

+ InoN2 + ^NoN3 + ^N2N3, (7) 

where the interaction energy is measured in units of vq. 
The diagrams shown in Fig. 2 lower the energy, as second- 
order perturbation theory implies. For example, for the 
process shown in Fig. 2(a) where two particles with m = 2 
scatter to a state with m = and m = 4, and then back 
to the initial state, the matrix element for each vertex is 
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M - ^^/N2iN2 - l)(iVo + l){Ni + 1) vo, 
Id 
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where Nm is the number of atoms with angular momen- 
tum mh. In Eq. (||) we have included the factor N4 + I for 
clarity, although as we mentioned above 7V4 is assumed 
to be equal to zero. The result of Eq. (^) is most easily 
derived by writing the interaction energy, Eq. (|2|) , as 
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where ai and a\ are annihilation and creation operators 
respectively, and 



i,j,k,l 



$*(r)$*(r)$fc(r)$((r)dr 
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when i+j ~ k-\-l and zero otherwise. In the case we con- 
sider here only the integrals /o.4,2,2 ^^nd Iifi.2,2 give a non- 
zero result reflecting the symmetrization of the bosonic 
wavef unction. 



As Eq. (^ implies, the difference between the energy 
in the intermediate and the initial states is IS.E = + 
Eq — 2e2 = — iVuo/8. Therefore, according to perturbation 
theory, the correction to the energy is 



AE 16 ^ ' ° 



(11) 



in the limit L <^ N. Table I gives the correction to the 
energy for the processes shown in Fig. 2. Adding these 
terms to 

£(0) 

, we see that the corrected interaction energy 
is, in units of vq and to leading order in N~^, 



£(1) = -No{No - 1) 
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To conserve particle number and angular momentum, we 
have the following constraints 



Nq + N2 + N3 = N and 2N2 + SN^ = L. 



(13) 



Using them, we express in terms of N,L, and N3, 
thus finding 

1 27 

= -iV(2iV-L-2) + — iV3(A^3-l). (14) 

The number 27/68 coincides with the numerical result 
1.59/4 reported in Ref. [^. Here can take all the 
non-negative integer values that are consistent with the 
constraints of Eq. ( [l3|) . Thus the number of excited states 
described by Eqs. (|13|) and ( p^ ) is equal to the integer 
part of L/6. These statements are exact asymptotically, 
i.e., for N 00, since there are processes which couple 
the quadrupolc m — 2 excitations with the octupole m — 
3 excitations and they contribute terms of order 1/A^ 
to the excitation energy. Concerning the yrast state, a 
consequence of Eq. ( |l4| ) is that its energy is given by the 
expression 



fo - -i2N-L-2), 
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in agreement with Refs. [|7 yiC| , |l2| |. Finally a result of 
Eq. (|lj) is that the yrast state for L <C iV is dominated 
by quadrupolc excitations, i.e., to leading order. 



N2 
N 



1 L 
2N' 



(16) 



In summary we have developed an effective theory 
which describes the ground state and the low-lying ex- 
cited states of a weakly-interacting Bose-Einstein con- 
densate under rotation, in the limit where the angu- 
lar momentum is much less than the number of atoms. 
This study demonstrates that there are low-lying excited 
states which differ by an energy of order vq (and not 
Nvq), and we have found agreement with a previous nu- 
merical study ll3] of the same problem. We should point 
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out that for all the values of L/N we have examined 
numerically, except the case L/N <C 1, we found that 
the low-lying excited states are separated from the yrast 
state by an energy of order iVuo, and in that respect the 
limit L <^ N seems to be unique. 
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FIG. 1. The six diagrams contributing to the bare inter- 
action energy £^^\ Eq. (^). The straight lines denote atoms 
with angular momentum given by the numbers written above 
the hnes; dashed lines denote the interaction. 
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Diagram 


Energy 


1(a) 


Ao(Ao-l)/2 


1(b) 


3A2(A2 - 1)/16 


1(c) 


5A3(A3-l)/32 


1(d) 


A0A2/2 


1(e) 


A0A3/4 


1(f) 


5A2A3/8 


2(a) 


-3A2(A2 - 1)/16 


2(b) 


-5A3(A3 - l)/544 


2(c) 


-A2A3/8 


2(d) 


-A2/2 


2(e) 


-3A3(A2 + l)/4 



TABLE L The contribution of the diagrams shown in 
Figs. 1 and 2 to the interaction energy (in units of vq). 



FIG. 2. The five (additional) diagrams contributing to the 
interaction energy Eq. (n3). 
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